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Abstract. 

We calculate the uvby light curves of the helium-strong chemically 
peculiar (CP) star HD 37776. Assuming that the chemical peculiar- 
ity influences the monochromatic radiative flux due to mainly the 
bound-free processes and using the model of the surface distribution 
of helium and silicon on HD 37776, as derived from spectroscopy, 
we calculate the emergent fluxes for a series of its rotational phases. 
We show that it is possible to consistently simulate light curves of a 
chemically peculiar star and that the basic properties of variability 
of helium CP stars can be understood in terms of the model of spots 
with a peculiar chemical composition. 
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1 Introduction 



The periodic light variations are a common feature among the magnetic CP stars. It is 
believed that most of these variations are caused by the rotational modulation of the observed 
radiative flux due to spots of a different peculiar chemical composition present on the stellar 



surface. In Krticka et al. (2004) we investigated the effect of the uneven surface distribution of 



helium on the monochromatic emergent flux and showed that the distribution of the flux may 
be changed by the bound-free processes. In this work we calculate the light curves of helium- 
strong CP star HD 37776 using the spectroscopically derived surface distribution of helium 



and silicon by Khokhlova et al. (2000) To study the influence of the chemical peculiarity 



on the observed light variability besides helium we selected silicon as: (i) it is overabundant 
for a large group of chemically peculiar stars that show light variations (e.g., Si, He- weak, 
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and He-strong stars), and these stars form a photometrically homogeneous group; (ii) a large 
overabundance of silicon is observed on the HD 37776 surface. 



2 Atmosphere models of HD 37776 

Given the map of the helium and silicon surface distribution we computed the emer- 
gent flux for the every surface element, sized 1° x 1°, and then integrated over the 
all visible surface. The fluxes were computed using SYNSPEC, the code for spectrum 
synthesis, for 8 x 10 LTE models computed with TLUSTY (|Hubeny fc Lanz 1995|) for 



different silicon and helium abundance from the grid [Si/H] = —3,-2,-1,0,1,2,2.5,3 
and [He/H] = — 5, — 4, — 3, — 2, — 1, 0, 1, 2, 2.5, 3 to cover the extent of values given by 



Khokhlova et al. (2000) The effective temperature T e g and surface gravity log g of the models 
are the same (see Tab. [1]). The atomic species and the number of their energy levels included 
in the model atmosphere calculations are listed in Tab. [2j 



Table 1: Parameters of HD 37776 





22 000 K 


Groote & Kaufmann(1982) 


\ogg 


4.0 


Groote & Kaufmann(1982) 


inclination % 


45° 




Khokhlova et al. (2000) 




spots parameters 






Khokhlova et al. (2000) 





Table 2: The atomic species and the number of their energy levels included in the model 
atmosphere calculations 



Ion 


Levels 


Ion 


Levels 


Ion 


Levels 


Ion 


Levels 


Hi 


9 


Nl 


21 


Nel 


15 


Sin 


40 


Hii 


1 


Nil 


26 


Nell 


15 


Sim 


30 


He I 


24 


Nm 


32 


Nem 


14 


Si iv 


23 


Hen 


20 


Niv 


1 


Ne IV 


1 


Siv 


1 


He ill 


1 


Oi 


12 


Mgl 


13 


Sn 


14 


Cl 


26 


On 


13 


Mgll 


14 


Sm 


10 


Cn 


14 


Om 


29 


Mgm 


14 


Siv 


15 


Cm 


12 


Oiv 


1 


Mg IV 


1 


Sv 


1 
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We calculate the flux from model atmosphere for given silicon and helium abundances 
H(X, [He/H], [Si/H]). The flux iJ c ([He/H],[Si/H]) in the color c is obtained as a convolution 
of H(X, [He/H], [Si/H]) with a Gauss profile peaked at the central wavelength of the color A c , 
with dispersion a c . The radiative flux in the color c from individual surface elements H C (Q) 
(Q are spherical coordinates on the stellar surface) is obtained by interpolation of fluxes from 
the model grid. Finally, the total observed radiative flux is calculated as the integral over all 
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Figure 1: The radiative fluxes from the HD 37776 atmosphere calculated with different 
silicon abundances and the observed averaged dereddened flux distribution according to 



Adelman & Pyper (1985) The radiative flux was smoothed to better demonstrate the light 
separation. In the uvby colors, the silicon spots are evidently brighter than the regions with 
low silicon abundance. Silicon abundances lower than the solar do not significantly change 
the radiative flux. 



the visible surface elements 

H c = J H c (tt)u(9) cos 6dtt (1) 

(6 is the angle between normal to the surface element and line of sight, and u(9) describes 
the adopted linear limb darkening) and the observed magnitude difference is then 

Am c = -2.5 log (J^J , (2) 

where H* et is a reference flux. This was performed for 36 rotational phases (some of them are 
shown on Fig. [2]) . 



3 The origin of light variations 

The plot of the radiative flux calculated for different values of the silicon abundance (Fig. [1]) 
shows that stellar surface in silicon rich regions is brighter in the uvby colors than in the silicon 
poor ones. The increase of the flux in the optical domain is caused by the flux redistribution 
from the far UV region mainly due to the enhanced silicon bound-free absorption. 



4 Elemental surface distribution 

At different rotational phases, surface elements with different silicon abundance are seen on 
the stellar disc. During the light maximum around the phase (j) = 0.75 (see also Fig. [3]) the Si 
rich regions are seen, which have large flux in the uvby (see Fig. [1]). During the light minimum 
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Figure 2: Variation of the silicon abundance on the surface of HD 37776 derived from spec- 



troscopy by Khokhlova et al. (2000) , smoothed for different rotational phases. The Si over- 
abundant regions are bright on the stellar surface (see Fig. [1]) and white on the graph; the 
Si underabundant regions are dark on the stellar surface and blue (dark) on the graph. 

around the phase <fi = 0.1 the Si poor regions are seen on the stellar disc, which have lower 
flux in the uv by. 

5 Predicted light variations 

Due to the uneven surface distribution of silicon (and partly also due to helium) the emergent 
monochromatic flux varies with the location on the visible disc. As shown in Fig. [TJ the Si 
rich regions are brighter in the uvby colors (the brightness temperature is higher), whereas 
Si poor regions are darker. 

Using the model atmosphere fluxes derived for the different helium and silicon abundances 
we calculated the theoretical light curves (see Fig. [3]) and compared them with the observed 
ones. The observed light variations can be explained as alternation of the emergent flux due 
to the uneven surface distribution of silicon and helium. 



6 Conclusions 

We are able to simulate the light variations of the He-strong CP star HD 37776 assuming 
that the light variations are due to the uneven surface distribution of silicon and helium. 
The silicon and helium spots modify the emergent flux mainly due to bound-free transitions. 
The predicted light curves reproduce the observed ones very well in their overall shape and 
amplitude. We stress that the theoretical light curves were obtained from the spectroscopical 
silicon surface distribution without using any free parameter. We were able for the first 
time (to our knowledge) to simulate the light curve of a CP star. Due to the satisfactorily 
good agreement between the observed and predicted light variations we conclude that this 
modelling is a very promising way towards the explanation of the light variations of chemically 
peculiar stars. 
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Figure 3: The predicted light variations of HD 37776 calculated using the silicon and he- 



lium surface distribution after Khokhlova et al. (2000) and compared with the observed ones 
(empty circles). 



operated at CDS, Strasbourg, France and On-line database of photometric observations ofmCP stars 
(jMikulasek et al. 2007ft . 
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